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Abstract. Polarimetric Synthetic Aperture Radar (PolSAR) images are 
establishing as an important source of information in remote sensing appli- 
cations. The most complete format this type of imaging produces consists 
of complex-valued Hermitian matrices in every image coordinate and, as 
such, their visualization is challenging. They also suffer from speckle noise 
which reduces the signal-to-noise ratio. Smoothing techniques have been 
proposed in the literature aiming at preserving different features and, 
analogously, projections from the cone of Hermitian positive matrices 
to different color representation spaces are used for enhancing certain 
characteristics. In this work we propose the use of stochastic distances 
between models that describe this type of data in a Nagao-Matsuyama- 
type of smoothing technique. The resulting images are shown to present 
good visualization properties (noise reduction with preservation of fine 
details) in all the considered visualization spaces. 
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1 Introduction 

Among the remote sensing technologies, PolSAR has achieved a prominent 
position. PolSAR imaging is a well-developed coherent and microwave remote 
sensing technique for providing large-scaled two-dimensional (2-D) high spatial 
resolution images of the Earths surface reflectivity; see Lee and Pottier [4] . 

The phenomenon speckle in SAR data hinders the interpretation these data 
and reduces the accuracy of segmentation, classification and analyses of objects 
contained within the image. Therefore, reducing the noise effect is an important 
task, and multilook processing is often used for this purpose in single-channel 
data. 

According to Lee and Pottier [3] , the principle to preserve the polarimetric 
signature and Polarimetric SAR image smoothing requires: (i) for each element of 
the image should be filtered in a way similar to multilook processing by averaging 
the covariance matrix of neighboring pixels; (ii) the filtering should be executed 
independently for each element of the covariance matrix; and (iii) homogeneous 
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regions in the neighborhood should be adaptively selected to preserve resolution, 
edges and the image quality. 

The statistical modeling provides a good support for the development of 
algorithms for interpreting PolSAR data efficiently, and for the simulation of 
plausible images. Frery et al. |2l3j introduce statistical tests for analyzing contrast 
in PolSAR images under the scaled multilook complex Wishart distribution, 
which has been successfully employed as a statistical model in such images for 
homogeneous regions. Frery et al. [3J derive several distances and tests for the 
complex Wishart model. 

This work presents a new smoothing process for PolSAR imagery based 
on stochastic distances and tests between distributions. This process, beyond 
reducing the noise effect, maintains geometric features of the PolSAR data. Vasile 
et al. [10] use a similar adaptive technique, but the decisions are based on the 
intensity information while we use the complete complex covariance matrix. 

The paper is organized as follows: In Section [2] we summarise the model 
for polarimetric data. Section [3] we describe the smoothing process for PolSAR 
images using stochastic distances between complex Wishart distributions, and 
the visualization of this kind of data. Results are presented in Section |4j while 
Section [5] concludes the paper. 

2 The Complex Wishart Distribution 

PolSAR imaging results in a complex scattering matrix, which includes intensity 
and relative phase data [3J. Such matrices have possibly four distinct complex 
elements, namely Syy, Svh, Shv, an d Shhi where H and V refer to the 
horizontal and vertical wave polarization states, respectively. The complex signal 
backscattered from each resolution cell is characterized by the p-tuple scattering 
matrix vector y, where p — 3 for a reciprocal medium (Svh = Shv)', see Ulaby 
and Elachi [9]. 

Thus, we have a scattering complex random vector y = [Svv, Svh, ShhY, 
where [•]* indicates vector transposition. In PolSAR data, the speckle might be 
modeled as a multiplicative independent zero-mean complex circular Gaussian 
process that modules the scene reflectivity |9l8j , whose probability density function 
is 

f(y,S) = n ^rexp{-y*U- 1 y}, 

where | • | is the determinant, the superscript V denotes the complex conjugate 
transpose of a vector, £ is the covariance matrix of y. The covariance matrix 
S, besides being Hermitian and positive definite, has all the information which 
characterizes the backscattering under analysis. 

Multilook processing is intended to enhance the signal-to-noise ratio, thus, is 
calculated the averaged over L ideally independent looks of the same scene. This 
results in the sample covariance matrix Z given by Z = L^ 1 Yle=i ViV*i, where 
L is the number of looks y^, for £={1,2,..., L}, and the superscript '*' denotes 
the complex conjugate transposition. 
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According to Anfinsen et al. [T] , Z follows a multilook scaled complex Wishart 
distribution, denoted by Z ~ W(£,L), Having £ and L as parameters, it is 
characterized by the following probability density function: 



r3L\ y/\L-3 

f z (Z';S,L)= I ;J T ^\^ exp{-£ tr^Z')}, 



where r 3 {L) = 7r 3 n?=o F(L — i), r(-) is the gamma function, tr(-) is the trace 



\z\ L r 3 (L) 
i), r(.) is 1 

operator, and the covariance matrix of Z is given by 
£ = E{yy*} 



(1) 



£{SiS!} £{SiS 2 *} £{SiS 3 *} 
E{S 2 St} E{S 2 S^} E{S 2 SI} 
E{S 3 Sl} E{S 3 S* 2 } E{S 3 S* 3 } 



where E{-} and the superscript '*' denote expectation and complex conjugation, 
respectively. 



3 Stochastic Distances Filter 

The proposed filter is based on stochastic distances and tests between distribu- 
tions 3 obtained from the class of (h, </>)-divergences. It employs in a modified 
Nagao-Matsuyama set of neighbors [5] , presented in Figure [l] 
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Fig. 1. Nagao-Matsuyama neighbourhoods. 



The filtering procedure consists in checking which regions can be considered 
as coming from the same distribution that produced the data in the central block. 
The sets which are not rejected are used to compute a local mean of covariance 
matrices. If all the sets are rejected, the filtered value is updated with the average 
on the central neighborhood around the filtered pixel. 

Each filtered pixel has a 5 x 5 neighborhood, within which nine overlapping 
areas are defined. Denote 0\ the estimated parameter in the central 3x3 neigh- 
borhood, and (#2, ■ ■ • , #<)) the estimated parameters in the eight remaining areas. 
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We estimate Oi — {Si) by maximum likelihood, assuming that the number 
looks L is known; details can be seen in Anfinsen et al. pQ. Based on a random 
sample of size n, let {Z l7 Z 2 , ■ ■ . , Z n }, the likelihood function related to the 
W(27, L) distribution is given by 

C ^ Z>) = { \Z\*rz{L) Y fll^l^expi-itr^-^')}. (2) 



Thus, the maximum likelihood estimator for S is £ = n _1 Si=i ^j- 

The proposal is based on the use of stochastic distances on small areas within 
the filtering window. Consider that Z\ and Z$ are random matrices defined on 
the same probability space, whose distributions are characterized by the densi- 
ties fz 1 (Z';0i) and /^(Z';^), respectively, where 0\ and Oi are parameters. 
Assuming that the densities have the same support given by the cone of Hermitian 
positive definite matrices A, the h-ip divergence between fz 1 and /z ; is given by 

D$ { z uZi ) = h ( fj{%$$) dtf). (3) 

where h: (0, oo) —> [0, oo) is a strictly increasing function with h(0) = and 
h'(x) > for every x £ R, and 0:(O, oo) — > [0, oo) is a convex function [7]. 
Choices of functions h and (f> result in several divergences. 

Divergences sometimes are not distances because they are not symmetric. A 
simple solution, described in [21316] . is to define a new measure df^ given by 

D'){Zx,Z i )+Dh(Z h Zi) 
d\{Z x ,Z i )= + K %> 2 ^ V l> . (4) 

Distances, in turn, can be conveniently scaled to present good statistical properties 
that make them suitable as test statistics 13 6: 



sl{e u e i ) = 2 ^d%{e l A), (5) 

where &\ and Oi are maximum likelihood estimators based on samples size m 
and n, respectively, and k = (h' (O)^")^ 1 . When B\ = Oi, under mild conditions 
6>^ l (<?i,0i) is asymptotically \m distributed, being M the dimension of 0%. Ob- 
serving S^(0i, Oi) — s, the null hypothesis 0\ — Oi can be rejected at level r] if 
P r (xlr > s ) — V- Details can be seen in the work by Salicru et al. [7]. 

Since we are using the same sample for eight tests, we modified the value of 
77 by a Bonferroni-like correction, namely, the Sidak correction, that is given by 
77 = f — (f — a) 1 /*, where t is the number of tests and, a the level of significance 
for the whole series of tests. 

Frery et al. [5] derived several distances for the W(S, L) model, the one 
presented in Equation ([!]) among them. The statistical test used in this paper 
was derived from the Hcllingcr distance, yielding: 
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PolSAR is used to measure the target's reflectivity with four polarization 
channel combinations (HH, HV , VH and VV), which can be expressed as a 
complex scattering matrix [4 9J. Transformations on these channels polarization 
makes it possible to visualize the PolSAR data as a color image. 

Two ways of visualizing the covariance matrix in false color are the Pauli (in 
the horizontal basis) and Sinclair decompositions. They assign \Sjjh ~ Svv\i 
\2Shv\ 2 and \Shh + SW| 2 , and |<SVv|, |25w| and |»Sfm| to the red, green and 
blue channels, respectively. 



4 Results 



The NASA/ Jet Propulsion Laboratory Airborne SAR (AIRSAR) of the San 
Francisco Bay was used for evaluating the quality of the procedure. The original 
polarimetric SAR data was generated with 4-looks and 900 x 1024 pixels. Figure [2] 
presents results in the Pauli decomposition, while Figure [3 shows their coun- 
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terparts using the Sinclair decomposition. Figures 2(a) and 
and |3(b)| respectively) show the original data set and a zoom. 

Figures [2(c) and 2(d) (Figures 3(c) and |3(d)| resp.) show the effect of the 
mean computed on windows of size 5x5 over the whole image. Albeit the noise 
reduction is evident, it is also clear that the blurring introduced eliminates useful 
information as, for instance, curvilinear details in the forested area. 

Figures 2(e) and 2(f) (Figures 3(e) and |3(f)| resp.) present the result of 
smoothing the origintal data set computing means which pass the Hellinger test 
at the level significance a = 80%. The noise effect is alleviated, c.f. the reduced 
graininess specially in the forest and over the urban areas, but fine details are 
more preserved than when the mean is employed. The directional selectiveness 
of the proposed filter retains linear structures as, for instance, the streets and 
the docks. Bright structures within the forest are enhanced, and their linear 
appearance is maintained. 



5 Conclusions 



The proposed technique combats the effect of speckle noise in all the areas of the 
test image. The filter is selective, producing stronger noise reduction in untextured 
areas, while it preserves fine details as linear structures and the forest texture. 
The proposal was compared with the simple mean filter using the decomposition 
process for PolSAR images. 

Next steps will be assessing quantitatively the proposal, using iterated filters 
(since the complex Wishart distribution is closed under convolutions), and esti- 
mating the equivalent number of looks in order to consider possible departures 
from the homogeneous model. 
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Fig. 2. PolSAR data on Pauli Decomposition. 



8 Torres, Medeiros & Frery 




Fig. 3. PolSAR data on Sinclair Decomposition. 



